Temperature dependent correlations in covalent insulators 
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Motivated by the peculiar behavior of FeSi and FeSb2 we study the effect of local electronic 
correlations on magnetic, transport and optical properties in a specific type of band insulator, 
namely a covalent insulator. Investigating a minimum model of covalent insulator within a single-site 
dynamical mean-field approximation we are able to obtain the crossover from low temperature non- 
magnetic insulator to high-temperature paramagnetic metal with parameters realistic for FeSi and 
FeSb2 systems. Our results show that the behavior of FeSi does not imply microscopic description 
in terms of Kondo insulator (periodic Anderson model) as can be often found in the literature, but 
in fact reflects generic properties of a broader class of materials. 
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The effect of local electronic correlations turning met- 
als into Mott insulators has been one of the major themes 
of condensed matter theory. Electronic correlations in 
band insulators have received much less attention and 
were largely limited to investigation of Kondo insulators 
(KI) described by the periodic Anderson model at half 
filling. The possibility of a transition from a band insu- 
lator (BI) to a Mott insulator (MI) in the ionic Hubbard 
model [lj attracted new attention to the correlated Bis. 
Among the correlated Bis FeSi 0], and more recently 
FeSb2 have been subject of a particular interest as 
they are the only known 3d materials exhibiting behav- 
ior similar to Kondo insulators. At low temperatures 
the electrical resistivity p(T) and magnetic susceptibility 
x(T) correspond to thermally activated narrow gap semi- 
conductor. However, at higher temperatures (T >300A") 
the slope of _p(T) changes sign and becomes positive as 
in a metal 4}. The metal-insulator crossover is consis- 
tent with the optical measurements where the gap in the 
low-T absorption spectrum starts getting filled with in- 
creasing temperature, with no trace of the gap observed 
0,0, E3] any more at T >300K. The magnetic suscepti- 
bility x(T) vanishes at low T, passes through a maximum 
near 500 K, and obeys a Curie- Weiss law at higher tem- 
peratures [3, H| ■ Low electron doping by substituting Fe 
by Co («10%), or the substitution of Si by Ge («30%) 
yield a ferromagnetic metal 0- Recently a colossal See- 
beck coefficient was reported in FeSb2 |8|. 

LDA band structure calculations predict FeSi to be 
a nonmagnetic semiconductor [9( in agreement with the 
low temperature experimental data. Interestingly, con- 
trary to the usual underestimation of the gap by LDA, 
the calculated value «100meV is 2-3 times larger than 
the estimate from optical and photoemission measure- 
ments 0, HI, El- Using the LDA+U approach to mimic 
the effect of on-site Coulomb repulsion, Anisimov et al. 
[irl | found a ferromagnetic metallic phase very close by in 
energy. This approach has also allowed to describe tran- 



sition from nonmagnetic insulator to ferromagnetic metal 
with substitution Ge for Si in isoelectron electron series 
FeSii--rG e,. Jill , and the semiconductor-metal transition 
in FeSb2 [12j |. Nevertheless, current bandstructure the- 
ory itself cannot explain the temperature dependence of 
physical properties. 

To describe the susceptibility and other thermody- 
manic quantities Mandrus et al. [13| proposed a temper- 
ature independent model density of states with two very 
narrow («500if ) peaks separated by a gap of (wlOOOif ). 
The origin of the narrow peaks was attributed to an 
extreme renormalization of the non-interacting bands 
similar to Ce-based Kondo insulators invoking localized 
states very weakly hybridizing with broad conduction 
band. However, the assumption of strongly renormal- 
ized peaks surviving up to high temperatures appears 
rather unlikely. Moreover, the Fe3d-states hybridize very 
strongly with sp-bands: a large contribution of Fe3d- 
orbitals character can be observed in rs8 eV energy region 

am. 

In this work we study a minimal model of correlated 
covalent insulator Jli| within the dynamical-mean field 
approximation (DMFT). While the construction of the 
model is guided by the bandstructure of FeSi and FeSb2 , 
we do not aspire to describe the physics of these mate- 
rials in its full width. Our goal is to use simple model 
free of the complexity of a multi-band low-symmetry sys- 
tem to capture semi-quantitatively the essential physics 
of FeSi/FeSb2 and provide a well understood reference 
point for later multi-band study. 

In construction of the model we envision Hubbard 
Hamiltonian on a lattice with one orbital per site, but 
more than one site in the primitive cell. The DMFT 
approximation allows us to calculate important phys- 
ical quantities, in particular the single-particle spec- 
trum, from the knowledge of the non-interacting spec- 
tral density D(uj) and the on-site repulsion U (U—1.5 
eV throughout this work). The key feature of D(w), ob- 
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FIG. 1: (color online) Single-particle spectral density at var- 
ious temperatures. The non-interacting spectral function is 
marked with the dashed-line. The inset shows the tempera- 
ture dependence of dc conductivity. 
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FIG. 2: (color online) Single-particle self-energy for var- 
ious temperatures. 



tained by truncating the LDA density of states for FeSb2 
so that the low energy peak and gap structure is cap- 
tured, is the presence of a gap between bands of the 
same orbital character. Fu and Doniach 15] studied a 
two-orbital model with an on-site hybridization, which 
leads to a gapped D{lo). Furmann et al. [16| calculated 
the single-particle and optical spectra of such model and 
addressed the question of interaction driven BI to MI 
transition at high temperature. 

In general, gapped local spectral function D(uj) is char- 
acteristic of covalent insulator, which is an instance of a 
band insulator. In contrast to ionic insulator, arising e.g. 
from staggered potential or crystal field, with the local 
orbitals either filled or empty, the local orbitals of cova- 
lent insulator are half (or in general partially) filled. The 
groundstate of a covalent insulator is characterized by 
formation of non-local bonding combinations of orbitals, 



which are separated by a gap from the anti-bonding com- 
binations. Introducing an on-site repulsion U leads to 
competition between electron localization and formation 
of non-local bonds. Large U makes the system a Mott 
insulator (MI) [17] . The concept of BI to MI transitions 
driven by the on-site repulsion was recently studied in the 
context of the ionic Hubbard model [l[ . While analogies 
can be drawn between the two cases there are also obvi- 
ous differences which make the details BI-MI transitions 
in the correlated covalent insulator an interesting ques- 
tion. Here we focus on evolution of physical properties 
with temperature. 

The central quantity in the DMFT framework is the 
single-particle self-energy which is determined by 

the self-consistent solution of the equations: 



G(iuj n ) = J dv 
G^\iu> n ) = G 



D(u) 



ILUn + (1 - V - S(iw n ) 

_1 (iw n ) - S(iw n ), 



(1) 
(2) 



where E[Go(zw n ), U] is determined by solving of the aux- 
iliary Anderson impurity problem, for which wc employ 
the Hirsch-Fye quantum Monte-Carlo algorithm [ijj]. In 
Fig. Q] we show the evolution of the single-particle spec- 
tral density with temperature T . At low T renormalized 
quasiparticle bands with reduced gap (relative to U = 
case) and Hubbard bands are well distinguished. With 
increasing temperature the lifetime of the quasiparticles 
decreases (see the self-energy in Fig. [2]) and the sharp 
peaks in the spectral density are smeared as observed 
in the photoemission experiment. [f| The system thus 
smoothly evolves from low-T insulator to high-T inco- 
herent metal. 

This behavior becomes apparent from the temperature 
dependence of dc conductivity <r(T) calculated according 
to 



a(T) oc / deduj<S>(e)p^(e,uj;T) 



df(u;T) 



p{e,u>\T) 



1 



S(«;T)-i)- 



dk x ' 



-D(e). 



(3) 
(4) 
(5) 



Approximating the squared electron velocity (ff^) 2 by 
its average value the transport density $(e) reduces to 
the spectral density D(e) times a constant. Equation 
(J3j) is only valid if the quasiparticle lifetime is finite. To 
ensure this at all temperatures we introduce a small con- 
stant contribution to the linewidth 1/r = 0.015 eV which 
can be thought of as arising from impurity scattering in 
real materials. Temperature enters formula ([3]) through 
the Fermi-Dirac distribution f(u>;T) and the self-energy 
H(u>;T). The T-dependence of resistivity p = er -1 (see 
Fig. QJ is governed by the Fermi-Dirac distribution at 
low temperatures where the thermally activated behav- 
ior of band insulator is observed. At higher temperatures 
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FIG. 3: (color online) Intra-band contribution to the optical 
conductivity at various temperatures. The inset shows the 
evolution of a typical inter-band contribution. 

the increase of the number of thermally activated charge 
carriers is compensated by the decreasing quasiparticlc 
lifetime, which eventually becomes dominant and the re- 
sistivity assumes the metallic T-dependence with positive 
slope. This behavior of resistivity is observed in both 
FeSi 0] and FeSb 2 @]. 

In general, the conductivity consists of intra-band and 
inter-band contribution. The intra-band contribution is 
obtained from [lj| 

cr(u)) oc / dedu)'Q(e)p(e, uj r )p(e, uj' + uj) 

(6) 

x ,/W)-/(^ + ^ 

UJ 

To calculate the inter-band contribution the band disper- 
sion is necessary. In the present case of a narrow band 
semiconductor the behavior of the inter-band contribu- 
tion can be understood by looking at a typical transition 
corresponding to a specific fc-point. We choose excitation 
from the valence band peak at t\ to conduction band 
peak at ti as the typical transition 

<jj B {uj,k) ocM k / duj' p(e 1 ,Lu')p(e2,u;' + uj) 

(7) 

/(iy)-iW + c) 

x . 

UJ 

The total IB conductivity <jib(w) would be obtained by 
summation over the k-points weighted by the transition 
amplitude M/-, which results in some broadening due to 
the band dispersion. The T-dependence of contributions 
from different fc-points follows the same pattern governed 
by the self-energy T) and the Fermi-Dirac function, 



and thus the typical transition provides a good idea about 
the behavior of total ctjb(w). The intra-band contribu- 
tion is shown in Fig. [3J The main effect of decreasing 
the temperature consists in depletion the low frequency 
region of spectral weight. Drudc peak appears at inter- 
mediate temperatures, due to the thermally populated 
quasiparticle states, and disappears at lower tempera- 
tures. The appearance of the Drude peak is not a generic 
feature of the model, but depends on the relative size 
of the Kondo scale, give by the interaction parameter 
U and the overall bandwidth, and the semiconducting 
gap. The appearance of well defined quasiparticle states 
at low temperatures is also reflected by building-up of 
a peak in <7i n ter (inset of Fig. [3|) above the energy of 
the single-particle gap. Absence of clear Drude peak at 
high temperatures, and spectral weight transfer from low 
frequency region to above the single-particle gap upon 
cooling agrees well with the trend observed experimen- 
tally 

Next we discuss the behavior of local and uniform mag- 
netic susceptibilities, which are calculated by integrating 
the imaginary time correlation function (m 2 (r)m z (0)) 
and employing the formalism of Ref. [20| respectively. In 
Fig. |4] susceptibilities of the system without local inter- 
action are compared will those of interacting and doped 
system. At half filling the uniform susceptibility exhibits 
an exponential decrease at low T, while the local sus- 
ceptibility saturates to a finite value, in both interacting 
and non-interacting system. The interaction, which fa- 
vors formation of local moments, strongly enhances both 
local and uniform susceptibility and results in appear- 
ance of Curie- Weiss tail at high T. Doping of 0.2 elec- 
tron has a profound effect on the susceptibility. At this 
doping level the chemical potential is located close to the 
valence band peak and thus a Stoner-like enhancement 
of the susceptibility can be expected. While the high T 
local moment is reduced due to departure from half fill- 
ing, at low T both the uniform and local susceptibility 
show an increase characteristic for local moment systems. 
Different slopes of the inverse local and uniform suscep- 
tibilities of the doped system (inset of Fig. [4|) indicate 
that a simple Heisenberg model with rigid local moments 
does not provide the appropriate picture. 

Taking the local dynamical correlations into account 
via DMFT method results in a satisfactory description of 
all major FeSi-FeSb2 anomalies. The energy gap present 
in the single-particle spectrum at low temperatures, re- 
duced by a factor of two from its LDA value, starts fill- 
ing with increasing T, accompanied by smearing of the 
sharp peaks, so that above 600K no trace of the gap is 
visible. Corresponding change from insulating to metal- 
lic character is observed in the resistivity and the op- 
tical conductivity. The uniform magnetic susceptibility 
increases exponentially from T = 0, reaches maximum at 
fz600K and approaches Curie- Weiss law at higher tem- 
peratures. We find a moderate (Z rj 0.5) very weakly T- 



4 




500 



1000 1500 
Temperature (K) 



2000 



FIG. 4: (color online) Temperature dependence of local 
(black) and uniform (red) susceptibility for half-filled (dia- 
monds) and doped (squares) system. The dashed lines mark 
the local (black) and uniform (red) susceptibility of the non- 
interacting system. 



dependent quasiparticle renormalization below the room 
temperature implying that in this temperature range the 
T-dependence of physical quantities such as conductivity 
or thermopower enters predominantly through the Fermi- 
Dirac distribution. This behavior is quite different from 
classical, e.g. Ce based, Kondo insulators and heavy- 
fermion compounds where a strong T-dependence of the 
self-energy and thus the quasiparticle renormalization ex- 
ists down to very low temperatures. Sharpening of the 
spectral features at low temperatures due to quasiparti- 
cle renormalization enhances the thermopower, propor- 
tional to the first moment of spectral density [19(, and 
thus qualitatively agrees with the experimental observa- 
tion y. 

In real materials with multiple bands the orbital occu- 
pations deviate from half-filling, e.g. due to crystal-fields. 
Therefore the effect of competition between covanlency 
and interaction-driven localization such as correlation en- 
hancement of susceptibility is expected to be smaller then 
in the studied model. One can also imagine a situation 
with a staggered potential and a site dependent U, which 
can then describe materials with strong metal-ligand hy- 
bridization such as LaCoOs. We believe that presented 
minimal model can provide a useful reference point for 
understanding such systems. 

In conclusion we have investigated a minimal model of 
correlated covalent insulator to study the unusual tem- 
perature evolution of physical properties of FeSi and 
FeSb2 . The model captures all the remarkable properties 
of these materials including the insulator to bad metal 
transition and the appearance of Curie- Weiss suscepti- 
bility at elevated temperatures. 
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